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Abstract. We study nonequilibrium properties of an atomic quantum 
dot (AQD) coupled to a Bose-Einstein condensate (BEC) within Keldysh- 
Green's function formalism when the AQD level is varied harmonically 
in time. Nonequilibrium features in the AQD energy absorption spec- 
trum are the side peaks that develop as an effect of photon absorption 
and emission. We show that atoms can be efficiently transferred from 
the BEC into the AQD for the parameter regime of current experiments 
with cold atoms. 



1 Introduction 

The tremendous experimental progress in the field of ultracold atoms has made possi- 
ble with unprecedent control the realization and design of many-body atomic physics 
by means of very stable optical lattices and by the use of Feshbach resonances [JJ [2]. 
That has triggered the study of such systems for quantum information processing 
purposes [3] and as quantum simulators of condensed matter Hamiltonians @]. In 
this context the experimental development of atom chips [5] opened new perspectives 
towards atomic mesoscopic physics. Indeed the possibility to build up atomic waveg- 
uides above microfabricated surfaces allows one to study a variety of phenomena such 
as the interplay between interaction and transport through point contacts [5] , the dy- 
namics of soliton-like structures in waveguides [7] and the coherent flow of a Bosc 
Einstein condensate (BEC) through a double barrier potential realized as a quantum 
dot in a magnetic waveguide [5]. 

Another route towards trapping and manipulation of single atoms [5] has been re- 
cently proposed, which relies on a focused laser beam superimposed to a trap holding 
an atomic Bose Einstein condensate; that leads to the formation of an atomic quan- 
tum dot (AQD) [10], i.e. a single atom in a tight trap, which is coupled to a superfluid 
reservoir via laser transitions. In particular, in the case of a one dimensional super- 
fluid reservoir, the system maps onto a spin-boson model with ohmic coupling, which 
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exhibits a dissipative quantum phase transition [TT] and provides a spectroscopic tool 
to engineer open quantum systems and study their irreversible dynamics |12j . Such 
a system can also be employed as a probe of BEC phase fluctuations via dephasing 
measurements of the internal states of the AQD, as discussed in [T3]. More recently 
an AQD coherently coupled by optical transitions to two BEC reservoirs has also 
been considered together with a Josephson tunneling between the two BECs [T3] . Fi- 
nally, a bosonic single-impurity Anderson model (AM) has been recently numerically 
investigated |1 5) in order to understand the local dynamics of an AQD coupled to 
a BEC system, accounting for the entanglement and the decoherence of the macro- 
scopic condensate. In general, the strong collisional interaction in the AQD can locally 
break a BEC state to bring up the excitations of normal particles inside the dot and 
thus a non-equilibrium particle transport can be realized, while the number of ex- 
cited particles in the QD gives an indirect probe of the coherence of the macroscopic 
condensate. 

In this work we discuss the realization of coherent particle transport in an AQD 
resonantly coupled to a BEC via a Raman transition with effective Rabi frequency 
and realized by harmonic time modulation of the AQD level. We do consider a fully 
quantum non-equilibrium situation within a Keldysh approach, not simply looking at 
the stationary limit, and focus on the local spectral properties of the QD. 



2 Model and results 

In our model system we consider cold bosonic atoms with two hyperfine ground states 
a and b. Atoms in state a, which form the reservoir, are confined in a shallow trapping 
potential V a (x) while atoms in state b are localized in a tightly confining potential 
Vb (x) . The atoms are identical bosons of mass m and form a Bose-Einstein condensate 
in a state a, i.e. a quantum reservoir. Furthermore states a and b are coupled via a 
Raman transition with effective Rabi frequency ft. In the collisional blockade limit 
for atoms b only a single atom is localized on the dot. Atoms a and b are described by 
the field operators W a and tf'j,, with 3^, (x) = ipb (x) b, where the operator b destroys 
an atom b in the dot in the lowest vibrational state with wave function tpb (x) and 

\P a ( x ) ~ Pa (x) 1 ^ 2 e~ l< ^ x ) in the long- wavelength limit, where p a and </> are the density 
and phase operator, respectively. The general many-body Hamiltonian of our system 
is [TD]: " 

H = H a + H b + H ab (1) 

where 



if, fh 2 



Ha= 2 dx {m ps 



V</> 



2 (x) + ^7T 2 (x) j . (2, 

Pa 



and 



2 



H b + H ab = (-Mo + g a b J dx | Vfo ( x ) | 2 Pa (x)^ b + b + 
b+b+bb + hQ{t) J dx cos(k ■ x) ($ a (x) §1 (x) e iw '* + h.c)j , (3) 



where the first term in Eq. ([3]) contains the Raman detuning 5q and the collisional 
interactions between the dot and the reservoir, g ab = ^ h r ^ ah , with effective mass 

m = m a m b /(m a + m b ) and s-wave scattering lengths a ab ] p a (x) = (x) W a (x) is 
the density operator for atoms a, which can be expressed in terms of the density 
fluctuation operator 77 (canonically conjugated to cf>) as p a (x) = p a + 77 (x); the 
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second term describes the interaction strength of two-body s-wave collision Ubb oc gbb 
between the dot atoms (which we do consider in the collisional regime) while the third 
one is the laser induced coupling between atoms a and b with Rabi frequency Q. We 
do consider a step-like pulse with frequency oj p and duration r. Since in the collisional 
blockade limit uj gap T ^> 1 where io gav is the gap between the single and two-atom, 
the dynamics describing the coupling to all other bound states can be neglected and 
the laser resonantly couples the condensate with the single atom ground states of 
the dot with an effective strength = TiQ J <ixcos(kx)-0 Q (x)-0£,(x), i.e. the Fourier 
transform of the wavefunctions overlap. At low enough temperatures, the reservoir 
Hamiltonian ([2]) is that of a Bose superfluid with equilibrium density p a , superfluid 
density p s and low energy phonon excitations of linear dispersion w = v s |q|, v s being 
the sound velocity. Eq. © can also be written as a bath of harmonic sound modes: 

H a = hv s £ |q| bp* (4) 
q 

where 6 q are standard phonon operators, defined by: 

?(x) = *£ q 
£(x) = £ q 

V being the system volume. The Hamiltonian (fTJ) can be reduced to that of a spin- 
boson model [TT] in the collisional blockade limit (i.e. large strength Ubb) when the 
internal state of the AQD is described by a pseudospin 1/2 with the spin-up and 
spin-down state corresponding to single and no atom occupation respectively. Thus 

the following replacements can be taken: b + b — > ^ 1+ 2^ an d b + — > <t + , where cr +tZ 
stand for spin operators. 

In order to describe the density fluctuations regime and to conserve the total 
number of bosons, the previous model can be conveniently mapped onto a bosonic 
single-impurity Anderson model (AM) [T3] : 

h = eb+b + jb+b (b+b - 1) + ek^&k + ^(*) ( 6+5k + b i b ) ■ 

k k 

where &k, b^ are annihilation and creation operators of noninteracting bosons confined 
in a shallow potential, e, U are immediately recognized as a function of the detuning 
and the collisional energies of ([5]). eu — v s \k\ and the last term is the laser- induced 
hybridization between particles in the AQD and the bosonic bath where Vk has been 
given above. An estimate of the Hamiltonian parameters is the following [1^]: the 
reservoir can be made of a condensate of N ~ 10 3 87 Rb atoms with density rib — 3T0 13 
atoms I cm 3 in a harmonic trap (with trapping frequency Vb = 2ir x 100-ffz), while the 
tweezers trap for the AQD in the collisional regime should have frequencies of the or- 
der of hundreds of kHz and we take v a = 2ir x 100kHz and a detuning 5q — 284kHz. 
The Raman coupling between the dot and the reservoir for these parameters can be 
estimated of the order of J? ~ 30kHz. 

Since we are interested in the nonequilibrium properties and particle transport 
through the AQD, we do consider a situation in which the AQD level is harmonically 
varied in time via the coupling to the radiation, i.e. e — > e(t) = eq + e^coswi, and 
employ the real-time Keldysh formalism to calculate the spectral properties of the 
AQD and particle transport through it. Since in the following we are interested in the 



2UqVp a 

2mv B V 



1/2 



2 iq X Oq - b- 

iq ' x (&q + bt 



(5) 
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collisional blockade limit of the AQD, the dot levels eo and have to be considered 
not as bare ones, but as renormalized by the requirements of a single atom occupation. 

The BEC energy absorption and the particle current flowing through the AQD 
are given by: 



E k (t)) = --([E k ,H(t)]), 



(?) 
(8) 



where — Eub^bk is the BEC Hamiltonian and N = b^b^. By computing 
the commutators, both the BEC energy absorption and the particle current can be 

fa- 



expressed in terms of the lesser Green's function G^ k (t,t ^ = —i (b^ (t'^j b (t) 



E k (t) ) = e k J k 



J it) 



k 



Jk 



where 



Jk 



-Re 



{Vk(t)G< M (t,t)} 



as: 



(9) 



(10) 



is the momentum resolved current. By means of the equations of motion method, a 
general expression for the contour ordered Green function can be found in terms of 
the AQD and BEC Green's function, G and <7k respectively: 



o,k (m') = f dfiG(Mi)tE(ti)flk(t 1 ,t / ) 



(11) 



By applying the Langreth rules |17) and using the explicit expression for the BEC 
Green's function, i.e. g< (t,t'^j = if be (ek) e"^'-''), and g£ a (t,A = 

= T i0 (±t =p *') e" iEk (*"*') , where f BE is the Bosc distribution function of the reser- 
voir, one can rewrite the momentum resolved particle current as: 



Mt) 



/m(y k (t) / dhV£ 

I J —oo 



iek(t-ti) 



[G r (i, t\ ) Jbe (£k) + G< (t,t x ) 



(12) 



In terms of the self-energies the momentum resolved particle current becomes: 



Mt) 



dti [G r (t,t 1 )S^ (t 1 ,t')+G < (t,ti)^k 



(13) 



where the advanced and lesser self-energies of the reservoir are (ti,t^j = 

= |Vk| 2 <?£ and £< (ti,t) =\V k \ 2 g< (*!,*'). 

When the dot level is varied adiabatically in time, i.e. e(t) = £o + e^ cos(wi) 
with <C Eo, the time dependence of the tunneling (t) can be neglected and 

where a and a are the ground state size of the dot and the 
BEC in the loosely confined direction, respectively. Let us note that, since a should 
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be larger than the mean free distance between the particles a 3> I and I ~ £ where 
£ is the healing length £ = h/mc, it follows that c > ( and we take it equal to 
10 _6 m while a is taken 10 _4 m t 13j. In general, for very tightly confined AQD the k 
dependence of Vk can be neglected. 

In this nonequilibrium problem all Green's functions and self-energies depend 
explicitly on two time variables, so one needs to resort to two time Fourier transform 
and Eq. (|T2|) becomes: 



Mt) = —Re | / dE ^E2dE 3 ei(E3 _ Ei)t [Cfr s< {E ^ Es) 

+G < (E 1 ,E 2 )ZZ(E 2 ,E 3 )]}, (14) 
The lesser and retarded Green function for the dot G <,r (E%, E 2 ) are given by: 

G<(E 1 ,E 2 )=J2 I ^^■G r {E 1 ,^)S<{^ 2 )G a ^,E 2 ), (15) 
k J ( 27T ) 

, 3 , 2 4 G r (E U E 3 ) SI (E 3 ,E 4 ) G r (E 4 , E 2 ) , 
where the zeroth-order dot Green's function in the collisional blockade limit is 

™ / i n , — i I dt e(t') —i I dt [eo+£oj cos(wt')l , ., 

G r (t,tx) = -i6(t-h)e Jt i v ' = -«0(t-ti)e J 'i L v n while the 

reservoir self-energies — ^ k E£ and S r = ^ k are, respectively: 

E< (E U E 2 ) = ^|Vk| 2 g< (E 1 ,E 2 ) 

k 

= -i (2tt) 2 J2 Jbe (£k) |14| 2 S (E 1 - e k ) S (e k - E 2 ) , (17) 

k 

S r (E 1 ,E 2 )=Y j m 2 9UEuE 2 ) 

k 

= - (2-) £ l^l 2 r {El ~+L = ir (Ei), (18) 

^ -E-i - £k + «0+ 

k 

where we have defined A(Ex) and r(Ei) as the principal value and imaginary part of 
l/(i?i — £k + i0 + ). By substituting the expression of G r and the self-energy in Eqs. 
(|15p and (|16|) and taking the limit << 1, and keeping only terms linear in w we 
get the following dot Green's functions, within single photon approximation: 

G r (Ex,E 2 ) = 2tt [AtS (Ex - E 2 ) + B+6 [Ex -E 2 +u) + B~5 (Ex -E 2 - u)] , 

(19) 

and 

G< (E X ,E 2 ) = ir (E 2 ) [f BE (Ex) AxA*8 (E 2 - Ex) + Ax (B+)* 6 (E 2 -E x +lj) + 

Ax (B 2 )*S (E 2 -Ei-u) + B+A^S (E 2 -Et-u) + 

Bi A* 2 S (E 2 - Ex + u)] , (20) 

where we have defined: 

a (M^)T o± "Ji(^)j°(^) 

^1,2 - Ei- S0 -A(Ei, 2 )+ir(Ei, 2 )> D l,2 - (E 1 ,2-e -A(E 1 , 2 )+'>.r(E 1}2 ))(E 1} 2±< J j-e -A(E 1}2 )+ l r(E 1}2 ))- 

(21) 
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Here Jo (77) and Ji (^ *fj are the 0-th and 1-th order Bessel functions, respectively. 
By inserting Eqs. (fT9]) - (l20"]) in Eq. (Tl4]) . we can numerically evaluate both the energy 
absorption rate (EAR) and the particle current. The behavior of the EAR and the 
particle current, within single photon approximation, are shown in Figs. 1 and 2. 

0.05 rr 

0.04; 

& 0.03 
< 

w 0.02; 
0.01 

0.00 L 

0.0 0.5 1.0 1.5 2.0 2.5 

Fig. 1. Energy absorption rate as a function of for the following values of parameters (in 
unit of the quantum dot energy level): eo = 1.0, ft = 0.1, e^ = 0.06, cj = 0.5, T = 1. The 
time is fixed at t = 1.0 in units of the inverse of uj. 




As shown, the EAR spectrum shows a peak at eo and two side peaks corresponding 
to the absorption and emission of a photon with energy uj. 
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The current induced by the quantum dot modulation follows the harmonic vari- 
ation of the perturbation and its amplitude and mean value get modified as a result 
of single-photon absorption and emission processes. 



3 Conclusions 

The nonequilibrium properties of an atomic quantum dot (AQD) coupled to a Bose- 
Einstein condensate (BEC) were studied within the Keldysh-Green's function formal- 
ism when the AQD level is varied harmonically in time. We have analyzed both the 
current and the AQD energy absorption spectrum in the full nonequilibrium situa- 
tion, within a single-photon approximation and shown that it is possible to achieve 
an efficient way of transferring particles from the BEC to the AQD in a range of 
parameters interesting for current experiments with cold atoms. Both quantities are 
actually a subject of active investigation theoretically as well as experimentally in 
order to probe and manipulate such systems. 

Indeed the coherent tunneling of particles from the BEC to the AQD could be 
employed to extract atoms on demand from the quantum reservoir, thus realizing a 
quantum tweezer [TB]. This high degree of control is needed in several protocols for 
quantum information processing with neutral atoms, in particular at the initialization 
stage of a quantum register |18j . Further applications can be envisaged, which run 
from the optimal control of atoms with microwave potentials for the implementation 
of quantum gates on an atom chip |19j to the realization of an efficient procedure to 
filter out from an optical lattice a preselected number of atoms per site |20) . In this 
way lattices with a desired site occupation could be engineered. 

On the other hand, the energy absorption rate could be employed as a spectro- 
scopic tool to probe several properties of the bosonic system under study, and in 
particular the current autocorrelation function, in analogy with a recent proposal 
by Giamarchi and co-workers [3T]. That allows one to directly probe the frequency 
dependent conductivity of the system. 
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